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ABSTRACT : Partial discharges (PDs) in high-voltage (HV) insulating systems originate from various local defects, 
which further results tit degradation of insulation and reduction in life span of equipment. One of the most widely used 
representations is phase-resolved PD (PRPD) patterns. For reliable operation of HV equipment, it is important to observe 
statistical characteristics of PDs and identify the properties of defect to ultimately determine the type of the defect. In this 
work, we have obtained and analysed combined use of PRPD patterns (f-q), (<p-n) and (n-q) using statistical parameters 
such as skewness and kurtosis for (f-q) and (<p-n),and mean, standard deviation, variance, skewness and kurtosisfor (n-q). 
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I. INTRODUCTION 

PD is an incomplete electrical discharge that occurs between insulation or insulation and a conductor. Partial 
discharges occur wherever the electrical field is higher than the breakdown field of an insulating medium. There are two 
necessary conditions for a partial discharge to occur in a cavity: first, presence of a starting electron to initiate an avalanche 
and second, the electrical field must be higher than the ionization inception field of the insulating medium [1]. In general, 
PDs are concerned with dielectric materials used, and partially bridging the electrodes between which the voltage is applied. 
The insulation may consist of solid, liquid, or gaseous materials, or any combination of them. PD is the main reason for the 
electrical ageing and insulation breakdown of high voltage electrical apparatus. Different sources of PD give different effect 
on insulation performance. The occurrence of sparks, arcs and electrical discharges is a sure indication that insulation 
problems exist. Therefore, PD classification is important in order to evaluate the harmfulness of the discharge [11]. 
PD classification aims at the recognition of discharges of unknown origin. For many years, the process was performed by 
investigating the pattern of the discharge using the well known ellipse on an oscilloscope screen, which was observed 
crudely by eye. Nowadays, there has been extensive published research to identify PD sources by using intelligent technique 
like artificial neural networks, fuzzy logic, and acoustic emission [11]. 

There seems to be an expectation that, with sufficiently sophisticated digital processing techniques, it should be 
possible not only to gain new insight into the physical and chemical basis of PD phenomena, but also to define PD 'patterns' 
that can be used for identifying the characteristics of the insulation 'defects' at which the observed PD occur [2]. 

Broadly, there are three different categories of PD pulse data patterns gathered from the digital PD detectors during 
the experiments. They are: phase-resolved data, time-resolved data and data having neither phase nor time information. The 
phase-resolved data consist of three-dimensional discharge epoch, cp charge transfer, q discharge rate, n patterns (<p~q, q~n 
and (p~n patterns) at some specific test voltage. The time-resolved data constitute the individual discharge pulse magnitudes 
over some interval of time, i.e., q~t data pattern. The third category of data consists of variations in discharge pulse 
magnitudes against the amplitude of the test voltage, V (for both increasing and decreasing levels), i.e., q~V data [3]. 

There are many types of patterns that can be used for PD source identification. If these differences can be presented 
in terms of statistical parameters, identification of the defect type from the observed PD pattern may be possible [4]. As each 
defect has its own particular degradation mechanism, it is important to know the correlation between discharge patterns and 
the kind of defect. Therefore, progress in the recognition of internal discharge and their correlation with the kind of defect is 
becoming increasingly important in the quality control in insulating systems [5]. Researches have been earned out in 
recognition of partial discharge sources using statistical techniques and neural network. In our study, we have tested various 
internal and external discharges like void, surface and corona using statistical parameters such as skewness and kurtosis for 
(cp-q) and (cp-n) and mean, standard deviation, variance, skewness and kurtosis for (n-q). 

II. STATISTICAL PARAMETERS 

The important parameters to characterize PDs are phase angle (p, PD charge magnitude q and PD number of pulses 
n. PD distribution patterns are composed of these three parameters. Statistical parameters are obtained for phase resolved 
patterns (cp-q), (cp-n) and (n-q). 

2.1 Processing of data 

The data to be processed obtained from generator includes tp, q, n and voltage V. From this data, phase resolved 
patterns are obtained. 
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2.1.1. Analysis of Phase-Resolved (<p-q) and (f-n) using Statistical Techniques 
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Fig.l (b) Block diagram of discharge analysis for (cp-n) 



PD pulses are grouped by their phase angle with respect to 50 (+ 5) Hz sine wave. Consequently, the voltage cycle 
is divided into phase windows representing the phase angle axis (0 to 360'). If the observations are made for several voltage 
cycles, the statistical distribution of individual PD events can be determined in each phase window. The mean values of 
these statistical distributions results in two dimensional patterns of the observed PD patterns throughout the whole phase 
angle axis [6]. A two-dimensional (2D) distribution cp-q and (p-n represents PD charge magnitude 'q' and PD number of 
pulses 'n' as a function of the phase angle '<p' [3]. 

The mean pulse height distribution H qn ((p) is the average PD charge magnitude in each window as a function of the 
phase angle (p. The pulse count distribution H n ((p) is the number of PD pulses in each window as a function of phase angle cp. 
These two quantity are further divided into two separate distributions of the negative and positive half cycle resulting in four 
different distributions to appear: for the positive half of the voltage cycle H qn + (cp) and H n + (cp) and for the negative half of the 
voltage cycle H qn " ((p) and H n " ((p) [5]. For a single defect, PD quantities can be described by the normal distribution. The 
distribution profiles of H qn (cp) and H n (cp) have been modeled by the moments of the normal distribution: skewness and 
kurtosis. 

ZlLiOq-lO'fOq) 



Skewness (S k ) = 



° 3 ZlLxfOq) 
gUxj ~ u) 4 f(x t ) 



■■(1) 



f(x) = PD charge magnitude q, 
u, = average mean value of q, 
o = variance of q. 

Skewness and Kurtosis are evaluated with respect to a reference normal distribution. Skewness is a measure of 
asymmetry or degree of tilt of the data with respect to normal distribution. If the distribution is symmetric, Sk=0; if it is 
asymmetric to the left, Sk>0; and if it is asymmetric to the right, Sk<0. Kurtosis is an indicator of sharpness of distribution. 
If the distribution has same sharpness as a normal distribution, then Ku=0. If it is sharper than normal, Ku>0, and if it is 
flatter, Ku<0 [3] [7]. 



2.1.2. Analysis of Phase-Resolved (q-n) using Statistical Techniques 
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Fig. 2 Block diagram of discharge analysis for (n-q) 
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Ku = kurtosis 

Statistical analysis is applied for the computation of several statistical operators. The definitions of most of these statistical 
operators are described below. The profile of all these discrete distribution functions can be put in a general function, i.e., 
y i= f(Xi). The statistical operators can be computed as follows: 

m u i t ^ X[Li(Xj) f ( x i) f~ 

Mean Value: (p) = (3) 

L i=\ i( x i) 

Variance: (o 2 ) = ^St^ (4) 

Standard Deviation = -^Variance (5) 

where, 

x = number of pulses n, 

f(x) = PD charge magnitude q, 

(i = average mean value of PD charge magnitude q, 

o = variance of PD charge magnitude q 

Skewness and Kurtosis are evaluated with respect to a reference normal distribution as described in section 2.1.1. 

III. RESULTS AND DISCUSSIONS 

Analysis involves determining unknown PD patterns by comparing those with known PD patterns such as void, 
surface and corona. The comparison is done with respect to their statistical parameters [9] [10]. 

3.1. Analysis for (<p-q) 

The phase resolved patterns are divided into two types: (cp-q) and (cp-n). The phase resolved patterns (<p-q) are 
obtained for three known PD patterns: void, surface and corona (as discussed in 3.1.1) and three unknown PD patterns: 
datal, data2 and data3 (as discussed in 3.1.3) [9] 



3.1.1. 2D distribution of ( D-q) for known PD patterns 

We have obtained the results from known PD parameters so as to plot the graphs showing below in Fig.3 (a), Fig. 3 
(b) and Fig.3 (c). These graphs are the phase cp vs. charge q plot for void, surface and corona discharges respectively. 




Fig.3 (a). Phase plot ((p-q) of void Fig.3(b).Phase plot ((p-q) of surface Fig.3(c). Phase plot ((p-q) of corona 

discharge discharge discharge 

3.1.2. Parameters of known PD patterns 

The table I below is the average value of skewness and kurtosis where values for H qn + ((p) are obtained over 0° to 
180° on the other hand values for H qn " (cp) are obtained over 181°to 360°, representing the phase. These values are obtained 
for known PD parameters. 



TABLE I. PARAMETERS OF KNOWN PD PATTERNS 



Parameter 


void 


surface 


corona 


Skewness H [|n + (cp) 


1.0013 


1.2134 


0.3555 


Skewness H qn " (cp) 


0.9901 


1.8219 


1.3659 


Kurtosis H qn + ((p) 


2.9046 


3.6064 


2.4354 


Kurtosis H m (cp) 


2.7872 


5.4506 


7.5947 
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3.1.3. 2D distribution of ( D-q) for unknown PD patterns 

We have obtained the results from unknown PD parameters so as to plot the graphs showing below in Fig.4 (a), 
Fig. 4 (b) and Fig.4 (c) which are the phase (p vs. charge q plot for datal, data2 and data3 respectively. 




Fig.4 (a) Phase plot (cp-q) of datal Fig.4 (b) Phase plot ((p-q) of data2 Fig.4 (c) Phase plot (cp-q) of data3 



From Fig.4 (a), it is seen that the following plot is similar to void and surface discharge. Fig.4 (b), is also similar to 
void and surface discharge and Fig.4(c), is similar to void discharge. 

3.1.4. Parameters of unknown PD Patterns 

The table II below is the average value of skewness and kurtosis where values for H qn + (cp) are obtained over 0° to 
180° on the other hand values for H qn " (cp) are obtained over 181°to 360°, representing the phase. These values are obtained 
for unknown PD parameters. 

TABLE II. PARAMETERS OF UNKNOWN PD PATTERNS 



Parameter 


datal 


data2 


data3 


Skewness H qn + (cp) 


0.8991 


0.7456 


1.0013 


Skewness H qll " (tp) 


1.1833 


0.8509 


0.9901 


Kurtosis H qn + (<p) 


2.5719 


2.1814 


2.9046 


Kurtosis H qn " (<p) 


3.7467 


2.6512 


2.7872 



3.2. Analysis for (<p-n) 

The phase resolved ((p-n) patterns consist of three known PD patterns: void, surface and corona (as discussed in 
3.2.1) and three unknown PD patterns: datal, data2 and data3 (as discussed in 3.2.3) [9]. The plots are discussed below: 

3.2.1 Phase resolved plot ( D-n) of known PD patterns 

We have obtained the results from known PD parameters so as to plot the graphs showing below Fig.5 (a), Fig.5 
(b) and Fig.5 (c) are the phase <p vs. number of pulses n for void, surface and corona discharges respectively. 




Fig.5(a).Phase plot ((p-n) of void Fig.5(b).Phase plot ((p-n) of surface Fig.5(c).Phase plot ((p-n) of corona 

discharge discharge discharge 

3.2.2. Parameters of known PD Patterns 

The table III below is the average value of skewness and kurtosis where values for H n + ((p) are obtained over 0° to 
180° on the other hand values for H n " ((p) are obtained over 181°to 360°, representing the phase. These values are obtained 
for known PD parameters. 
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TABLE III. PARAMETERS OF KNOWN PD PATTERNS 



Parameter 


void 


surface 


corona 


Skewness H n + (cp) 


0.4954 


1.0082 


1.3942 


Skewness H„" ((p) 


0.4329 


2.3686 


1.3798 


Kurtosis H„ + (cp) 


2.0535 


2.871 


4.8337 


Kurtosis H n " (cp) 


1.9137 


8.4788 


7.3215 



3.2.3. Phase resolved plot ( D-n) of unknown PD patterns 

We have obtained the results from unknown PD parameters so as to plot the graphs showing below Fig. 6(a), 
Fig.6(b) and Fig. 6 (c) are the phase (p vs. number of pulses n for void, surface and corona discharges respectively. 




Fig.6(a).Phase plot (cp-n) of datal Fig.6(b).Phase plot (cp-n) of data2 Fig.6(c).Phase plot (<p-n) of data3 



From Fig. 6 (a), it is seen that the following plot is similar to void and surface discharge. Fig.6 (b), is similar to void 
discharge and Fig.6 (c), is also similar to void discharge 

3.2.4. Parameters of unknown PD patterns 

The table IV below is the average value of skewness and kurtosis where values for H n + (cp) are obtained over 0° to 
180° on the other hand values for H n " (cp) are obtained over 181°to 360°, representing the phase. These values are obtained 
for unknown PD parameters. 

TABLE IV. PARAMETERS OF UNKNOWN PD PATTERNS 



Parameter 


datal 


data2 


data3 


Skewness H„ + (cp) 


0.8016 


0.574 


0.4954 


Skewness H„~ (cp) 


1.0169 


0.42 


0.4329 


Kurtosis H n + (cp) 


2.3724 


2.1091 


2.0535 


Kurtosis H„" ((p) 


3.2011 


1.8003 


1.9137 



3.3. Analysis for (n-q) 

The phase resolved patterns n-q are obtained for three known PD patterns: void, surface and corona (as discussed in 3.3.1) 
and three unknown PD patterns: datal, data2 and data3 (as discussed in 3.3.3) [10]. 

3.3.1. 2D distribution of n-q for known PD patterns 

We have obtained the results from known PD parameters so as to plot the graphs showing below Fig. 7(a), Fig. 
7(b), Fig. 7(c), Fig. 7(d) and Fig. 7(e) are the n-q plot of mean, standard deviation, variance, skewness and kurtosis for void 
discharge respectively. 




JtL. 



Fig.7(a).Mean plot (n-q) of void 
discharge 



Fig.7(b).Standard deviation plot 
void discharge 
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Fig.7 (d) Skewness plot (n-q) of void discharge Fig.7 (e) Kurtosis plot (n-q) of void discharge 



Referring to Fig. 7 (a), Fig. 7 (b) and Fig. 7 (c) of void discharge, it can be seen there is a peak occurring 
somewhere after 1500 cycle, which is a void discharge and in Fig. 7 (d) and Fig. 7 (e) of skewness and kurtosis, the value 
decreases at that cycle where peak occurs. 



We have obtained the results from unknown PD parameters so as to plot the graphs showing below Fig. 8(a), Fig. 
8(b), Fig. 8(c), Fig. 8(d) and Fig. 8(e) are the n-q plot of mean, standard deviation, variance, skewness and kurtosis for 
surface discharge respectively. 




Cycles 00 Cycles [n) 

Fig.8 (d) Skewness plot (n-q) of surface Fig.8 (e) Kurtosis plot (n-q) of surface 

discharge discharge 



In surface discharge, charges are distributed uniformly over all cycles for mean, standard deviation, variance, 
skewness and kurtosis as shown in Fig. 8(a), Fig. 8(b), Fig. 8(c), Fig. 8(d) and Fig. 8(e). Fig. 9(a), Fig. 9(b), Fig. 9(c), Fig. 
9(d) and Fig. 9(e) are the n-q plot of mean, standard deviation, variance, skewness and kurtosis for corona discharge 
respectively. 
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Fig.9 (a) Mean plot (n-q) of corona 
discharge 



Fig.9 (b) Standard deviation (n-q) of Fig.9 (c) Variance plot (n-q) of 
corona discharge corona discharge 




Fig.9 (d) Skewness plot (n-q) of corona 
discharge 



Fig.9 (e) Kurtosis plot (n-q) of corona 
discharge 



Referring to Fig. 9(a), Fig. 9(b) and Fig. 9(c) of corona discharge, it can be seen the charges starts occurring after 
500 cycle increasing somewhere upto 1200 cycle and then decreasing after 2000 cycle, and in Fig. 9(d) and Fig. 9(e) of 
skewness and kurtosis, the value decreases from 500 cycle till 2000 cycle. 

3.3.2. Parameters of known PD patterns 

The table V below is the average value of mean, standard deviation, variance, skewness and kurtosis for void, 
surface and corona discharges respectively . These values are obtained for known PD parameters. 

TABLE V. PARAMETERS OF KNOWN PD PATTERNS 



Mean 

Standard devia 



13320.32 
7553.716 



3.66*10" 17 



145.706 
126.009 



1.426 
1.139 



3.3.3. 2D distribution of (n-q) for unknown PD patterns 




Fig.10 (a) Mean plot (n-q) of datal 



Fig.10 (b) Standard deviation plot (n-q) 
of datal 



Fig.10 (c) Variance plot (n-q) of 
datal 
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Fig.10 (d) Skewness plot (n-q) of datal Fig.10 (e) Kurtosis plot (n-q) of datal 



We have obtained the results from known PD parameters so as to plot the graphs shown above Fig. 10(a), Fig. 
10(b), Fig. 10(c), Fig. 10(d) and Fig. 10(e) are the n-q plot of mean, standard deviation, variance, skewness and kurtosis for 
datal respectively. 

In Fig. 10(a), Fig. 10(b), Fig. 10(c), Fig. 10(d) and Fig. 10(e), the charges are uniformly distributed similar to 
surface discharge. Hence, it can be concluded that datal is having surface discharge. Fig. 11(a), Fig. 11(b), Fig. 11(c), Fig. 
1 1(d) and Fig. 1 1(e) are the n-q plot of mean, standard deviation, variance, skewness and kurtosis for data2 respectively. 




Fig.ll (d) Skewness plot (n-q) of data2 Fig.ll (e) Kurtosis plot (n-q) of data2 



In Fig. 11(a), Fig. 11(b), Fig. 11(c), Fig. 11(d) and Fig. 11(e), the charges are uniformly distributed similar to 
surface discharge. Hence, it can be concluded that data2 is having surface discharge. 




Fig.12 (a) Mean plot (n-q) of data3 Fig.12 (b) Standard deviation plot Fig.12 (c) Variance plot (n-q) of 

(n-q) of data3 data3 
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Fig.12 (d) Skewness plot (n-q) of data3 Fig.l2(e).Kurtosis plot (n-q) of data3 



Fig. 12(a), Fig. 12(b), Fig. 12(c), Fig. 12(d) and Fig. 12(e) are the n-q plot of mean, standard deviation, variance, 
skewness and kurtosis for data3 respectively. In Fig. 12(a), Fig. 12(b) and Fig. 12(c), there is a occurrence of peak after 1500 
cycle and in Fig. 12(d) and Fig. 12(e), the skewness and kurtosis value decreases at that peak which is similar to void 
discharge. Hence, it can be concluded that data3 is void discharge. 

3.3.4. Parameters of unknown PD patterns 

The table VI below is the average value of mean, standard deviation, variance, skewness and kurtosis for void, 
surface and corona discharges respectively . These values are obtained for unknown PD parameters 



TABLE VI. PARAMETERS OF UNKNOWN PD PATTERNS 



105.119 553.93 



16714.23 4.94*10 5 



1.64* 10 8 



IV. OBSERVATIONS 

From the above results, following observations are made: 

• Fig. 13(a), Fig. 13(b) and Fig. 13(c) are the characteristics of skewness and kurtosis (H qn + (cp) and H qn ~(<p)) of data 1, 
data2 and data3 against void., surface and corona discharges respectively. 

• From Fig. 13(a), it is observed that data3 characteristics overlaps void discharge characteristics, it can be concluded 
that data3 is void discharge. Data2 characteristics approximately fits against void, it can be concluded that data2 is 
also void discharge 

• From Fig. 13(b), it is observed that data 1 characteristics is similar to surface discharge characteristics, it can be 
concluded that data 1 is surface discharge. 

• From Fig. 13(c), it is observed that none of the data characteristics is similar to corona discharge characteristics, it 
can be concluded that none of the data has corona discharge. 




1 2 3 4 



Fig.13 (a) Characteristics of skewness 
and kurtosis(H qn + ((p) and H qn "(cp)) of 
datal, data2, and data3 against void 




12 3 4 



Fig.13 (b) Characteristics of 
skewness and kurtosis(H qn + ((p) and 
Hqn (cp)) of datal, data2, and data3 
against surface 



s 



4 




12 3 4 



Fig.13 (c) Characteristics of skewness 
and kurtosis(H qll + ((p) and H qn "(cp)) of 
datal, data2, and data3 against corona 
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Fig. 14(a), Fig. 14(b) and Fig. 14(c) are the characteristics of skewness and kurtosis (H n + ((p) and H n "(cp)) of data 1, 
data2 and data3 against void., surface and corona discharges respectively. These results are obtained over 0° to 180° for 
H n + (cp ) on the other hand values for H n " (cp) are obtained over 181°to 360°, representing the phase. These values are 
obtained for unknown PD parameters. 




Fig.14 (a) Characteristics of skewness Fig.14 (b) Characteristics of Fig.14 (c) Characteristics of skewness 

and kurtosis(H n + (cp) and H n "(cp)) of datal, skewness and kurtosis(H n + ((p) and and kurtosis(H n + ((p) and H n "((p)) of datal, 
data2, and data3 against void H„"(cp)) of datal, data2, and data3 data2, and data3 against corona 

against surface 



From Fig. 14(a), it is observed that data3 characteristics overlaps void discharge characteristics, it can be concluded 
that data3 is void discharge. Data2 characteristics approximately fits against void, it can be concluded that data2 is also void 
discharge. From Fig. 14(b), it is observed that datal characteristics are close to surface discharge characteristics, it can be 
concluded that data 1 is surface discharge. From Fig. 14(c), it is observed that none of the data characteristics is similar to 
corona discharge characteristics, it can be concluded that none of the data has corona discharge. 



In below figures, Fig. 15(a) represents the statistical characteristics of mean, standard deviation, variance, skewness 
and kurtosis of void discharge against data3. Fig. 15(b), Fig. 15(c) are the statistical characteristics of mean, standard 
deviation, variance, skewness and kurtosis of surface discharge against datal and data2 respectively. 




-ivvvvvv 12 5 4 5 12345 12345 



Fig.15 (a) Statistical Characteristics of Fig.15 (b) Statistical characteristics of Fig.15 (c) Statistical characteristics of 
data3 against void discharge datal against surface discharge data2 against surface discharge 

Plotting statistical parameters of void discharge against data3 in Fig. 15(a) shows data3 characteristics overlaps void 
characteristics, it can be concluded that data3 is void discharge. Similarly, for surface discharge, datal and data2 
characteristics (Fig. 15(b) and Fig. 15(c)) approximately fit surface discharge characteristics, it can be concluded that datal 
and data2 is surface discharge 

V. CONCLUSION 

From all the above observations using all the three phase-resolved patterns(cp-q), (<p-n) and (n-q), it can be 
concluded that, datal is surface discharge in all the three phase-resolved patterns, data2 is surface discharge in (n-q) pattern 
and void discharge in (cp-q) and (cp-n), hence data2 is both void and surface discharge and data3 is void discharge in all the 
three phase-resolved patterns. 
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